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ABSTRACT 

Context. Thermohaline mixing has recently been proposed to occur in low-mass red giants, with large consequence for the chemical 
yields of low-mass stars. 

Aims. We investigate the role of thermohaline mixing during the evolution of stars between 1 Mq and 3 Mq, in comparison with other 
mixing processes acting in these stars. 

Methods. We use a stellar evolution code which includes rotational mixing, internal magnetic fields and thermohaline mixing. 
Results. We confirm that during the red giant stage, thermohaline mixing has the potential to decrease the abundance of ^He, which 
is produced earlier on the main sequence. In our models we find that this process is working on the RGB only in stars with initial 
mass M < 1.5 Mq. Moreover we report that thermohaline mixing is also present during core He-burning and beyond, and has the 
potential to change the surface abundances of AGB stars. While we find rotational and magnetic mixing to be negligible compared to 
the thermohaline mixing in the relevant layers, the interaction of thermohaline motions with the differential rotation may be essential 
to establish the timescale of thermohaline mixing in red giants. 

Conclusions. To explain the surface abundances observed at the bump in the luminosity function, the speed of the mixing process 
needs to be more than two orders of magnitude higher than in our models. However it is not clear if thermohaline mixing is the only 
physical process responsible for these surface-abundance anomalies. Therefore it is not possible at this stage to calibrate the efficiency 
of thermohaline mixing against the observations. 

Key words. Stars: low-mass ~ Stars: evolution - Stars: abundances - Stars: interiors - Stars: magnetic field - Stars: rotation 



1. Introduction 

Stars are rotating, self-gravitating balls of hot plasma. Due 
to thermonuclear reactions in the deep stellar interior, stars, 
which presumably start out chemically homogeneous, develop 
chemical inhomogeneities. At the densities typically achieved 
in stars, thermal diffusion, or Brownian motion, is not able to 
lead to chemical mixing. However, various turbulent mixing pro- 
cesses are thought to act inside stars, leading to transport of 
chemical species, heat, angular momentum, and magnetic fields 
(lPinsonneaulti[T997t iHegeretaDIIOOOl I200I . 

Thermohaline mixing is usually not considered to be an 
important mixing process in single stars, because the ashes 
of thermonuclear fusion consist of heavier nuclei than its 
fuel, and stars usually burn from the inside out. The con- 
dition for thermohaline mixing, however, is that the mean 
molecular weight (jj) decreases inward. This can occur in 
accreting binaries, and the importance of thermohaline mix- 
' ing has long been recognized by the b i nary community (e.g ., 
Ide Greve (fe Cugierl [198 9: Sarn al 119921: IWellstein et aP i2001h . 



Recendv lCharbonnel & ZahnI (12007b'. CZ07) identified thermo- 
haline mixing as an important mixing process, which signifi- 
cantly modifies the surface composition of red giants after the 
first dredge-up. The work by CZ07 was initiated by the paper 
of Hggleton et al. (2006, EDL06), who found an mean molecu- 
lar weight (fi) inversion — i.e., ( p ) < — below the red 
giant convective envelope in a ID-stellar evolution calculation. 
While EDL06 then investigated the stability of the zone contain- 
ing the yu-inversion with a 3D hydro-code and found these layers 
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to be Rayleigh-Taylor-unstable, CZ07 could not confirm this, but 
found the layers to be unstable due to thermohaline mixing. 

lEggleton et al.l (l2006h found a //-inversion in their 1 Mq stel- 
lar evolution model, occurring after the so-called luminosity 
bump on the red giant branch, which is produced after the first 
dredge-up, when the H-burning shell source enters the chemi- 
cally homogeneous part of the envelope. The //-inversion is pro- 
duced by the reaction ^He(^He,2p)'*He, as predicted by UlricB 
( 119721) . It does not occur earlier, because the magnitude of the 
//-inversion is small and negligible compared to a stabilizing //- 
stratification. 

Mixing processes below the convective envelope in models 
of low-mass stars turn out to be essential for the prediction of 
their chemical yield of ^He (EDL06), and are essential to un- 
derstand the surface abundances of red giants — in particular 
the '^C/'^C ratio, ^Li and the carbon and nitrogen abundances 
(CZ07). This may also be important for other occurrences of 
thermohaline mixing in stars, i.e., in single stars when a //- 
invers ion is produ ced by off-center ignition in semi-degenerate 
cores (ISiessl l2009') or in stars which accrete chemically enriched 
matter from a companion in a close binary (e.g., Stancliffe et al] 
2007). Accreted metal-rich matter during the phases of planetary 
formation also leads to thermohaline mixing, which can recon- 
cile the observed metallicity distr ibution of the central stars of 
planetary systems jVauclaiHl2004{) . 

In the present paper we investigate the evolution of solar 
metallicity stars between 1 and 3 Mq from the ZAMS up to 
the thermally-pulsing AGB stage, based on models computed 
during the last years. We show for which initial mass range and 
during which evolutionary phase thermohaline mixing occurs 
and what consequences it has. Besides thermohaline mixing, our 
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models include convection, rotation-induced mixing, and inter- 
nal magnetic fields, and we compare the significance of these 
processes in relation to the thermohahne mixing. 



2. The speed of thermohaline mixing 

2.1. Thermohaline mixing VS Rayleigh-Taylor 

As pointed out by CZ07, a /i-inversion inside a star should give 
rise to thermohaline mixing, which is a slow mixing process act- 
ing on the local thermal timescale. Could a Rayleigh-Taylor in- 
stability be present in these layers? Indeed, EDL06 interprete 
the origin of the instability which they find in their 3D models 
as due to the buoyancy g produced by the /^-inversion, i.e. a 
dynamical effect. But a dynamical instability should only occur 
if the //-inversion were to lead to a density inversion. However, 
this would only be possible if the considered layers were convec- 
tively unstable in the hydrostatic ID stellar evolution models. 
As pointed out by CZ07 and as confirmed by our models, the 
yu-inversion produced by the ^He(^He,2p)'*He reaction does not 
induce convection. We conclude that the Rayleigh-Taylor insta- 
bihty may not be a likely explanation of the hydrodynamic mo- 
tions found by EDL06. A similar con clusion was also reached by 
iDenissenkov & Pinsonneaiil^ (l2008h . who studied in detail the 
instability driven by the //-inversion in both the adiabatic and the 
radiative limit. 

Because smaller blo bs have a smaller thermal timescale 
dKippenhahn et alj Il980l) . could EDL06 have found the high 
wavenumber tale of the thermohaline instability? They found the 
instability to occur within 2000 s. The size of a blob with such 
a short thermal timescale above the H-burning shell of a red gi- 
ant is on the order of 50 km. This is too small to be resolved in 
the 3D model shown by ELD06. Furthermore, an inspection of 
their Fig. 5 reveals that the length scale of the instability they 
found is on the order of 1 0^ ... 1 0'* km, which corresponds to ther- 
mal timescales of about 1 yr Therefore, it seems unlikely that 
EDL06 actually picked up the thermohaline instability in their 
3D hydrodynamic model, unless its non-linear manifestation in- 
volves a timescale much shorter than the thermal timescale. 

2.2. The efficiency of thermohaline mixing 

In Sect. |3] we explain the details of our implementation of ther- 
mohaline mixing in ID stellar evolution calculations. The dif- 
fusion coefficient for the mixing process contains a parameter 
which depends on the geometrical configuration of the fluid ele- 
ments. This parameter (ofth) is very important to understand the 
role of thermohaline mixing. It determines the timescale of the 
mixing (the velocity of the fingers/blobs) that we show in Sect.|4] 
and |5] plays a role not only in determining how fast the surface 
abundances of redgiants can change, but also if thermohaline 
mixing is present in stars of different mass and at different evo- 

luti onary phases. 

ICharbonnel & ZahnI (l2007bh showed that a high value of ath 
is needed to match the surface abundances of field stars after 
the luminosit y bum p (Grattonet al. 2000). Similar to the value 
adopted by lUMchl (i 19721) they use an efficiency factor corre- 
sponding to ath - 667 in our diffusion coefficient. This value 
corresponds to the diffusion proc ess involving fingers with an 
aspect ratio (length/width) of 5. 'Denissenkov & Pinsonneaultl 
(l2008h claim that to explain the observed mixing pattern in low- 
mass RGB stars, fluid elements have to travel over length scales 
exceeding their diameters by a factor of 10 or more. 



On the other hand, the order of magnitude of our effi- 
ciency parameter Cth = 2.0 corresponds to the prescription of 
Kippenhahn et al. (1980), where the diffusion process involves 
blobs of size L traveling a mean free p ath L before dissolving . 
The same prescription has been used bv lStancliffe et alj (l2007h . 
who dealt with the problem of thermohaline mixing in accreting 
binaries. The sensitivity of thermohaline mixing to a change of 
the efficiency parameter is shown in Fig.[T] where the change of 
the surface abundance of ^He at the luminosity bump is shown 
for different values of ath. 

That a prescription can reproduce the observed surface abun- 
dances may not be sufficient to prefer it over others. It is pos- 
sible that other mixing processes are at work. The resulting 
observed abundances could still be mainly due to thermoha- 
line mixing, as proposed by CZ07, but at t his stage it is not 
possible to exc lude that magnetic buoyancy dBusso et alj|2007l : 
Nordha us et alj |2008), or the interaction of different mixing 
processes (e.g., thermoh aline mixing and magnetic buoyancy, 
IDenissenkov et a n i2009l) . could play an even more important 
role. 

To clarify the picture here we discuss the main differences 
between the two physical prescriptions for thermohaline mix- 
ing. Experiments of thermohaline mixing show slender fingers 
in the lin ear regi r ne (e.g ., Krishnamurti 2003), supporting the 
picture of lUlrichI d 19721) . On the other hand the physical con- 
ditions inside a star are quite different from those in the lab- 
oratory. In particular the Prandtl number cr, defined as the ra- 
tio of the kinematic viscosity v to the thermal diffusivity /cj-, is 
very small in stars (cr ~ 10"''). This number is about 7 in water, 
where most of thermohaline mixing experiments have been per- 
formed. The question arises if for such small values of cr a finger- 
like structure can be stable, especially in layers where shear and 
horizontal turbulence is present. 2D hydrodynamic simulations 
of double-diffusive phenomena, with out external perturbations 
have b een performed in the past (e.g.. lMerrvfieldl[T995l : iBascoull 
l2007allbl) . While the resolution required by the physical condi- 
tions in stellar interiors is computationally not accessible, lower- 
ing the Prandtl number to values of abou t 10~^ always results 
in incr easingly unstable structures (Merr vfieldlll995l; IBascoull 
l2007ah . Therefore it might be dangerous to assume that the same 
configuration of thermohaline mixing, as observed in water, is 
occurring in stars. 

As we show in Sect. |6l the radiative buffer between the 
H-burning shell and the convective envelope is a region of 
the star where the angular velocity is rapidly changing. Even 
if the order of magnitude of the rotationally-induced instabil- 
ities is much lower than the one from thermohaline mixing 
(cf. Fig |9]), it is possible that the interaction of the shear mo- 
tions with the thermohaline diffusion prevents a relatively or- 
dered flow to be stable, contrary to Ulrich's assumption. That 
shear can decrease the efficiency of th ermohaline mixing was 
already pointed out by ICanutol (Il999h . The effect of strong 
horizontal tur bulence in stellar l a yers h as also been discussed 
by Denissen kov & Pinsonneaul3 (l2008l) . and we conclude in 
Sect. 16.31 that this effect could work against the fingers and in 
favor of blobs. 

The thermohaline mixing prescription proposed by 
^Kippenhahn e t al.l (Il980l) has been used for the calculations 
presented here. It is clear from the discussion above that a better 
calibration of the mixing speed requires realistic hydrodynamic 
calculations of the instability. 
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Fig. 1. Evolution of the surface abundance of ^He for o-,/, - 
(blue solid line), a,i, - 2 (red dotted line) and or,/, - 200 (green 
dashed line) from before the onset of thermohaline mixing to the 
core He-flash in a 1.0 Mq star. 

3. Method 



We u se a 1-D hydrodynamic stellar evolution code dYoon et alj 
I2OO6I and references therein). Mixing is treated as a diff'usive 
process and is implemented by solving the diffusion equation 



( dX„ \ 
I dt I 



d_ 

dm 



\ dm 



dt 



(1) 



where D is the diff'usion coefficient constructed from the sum of 
individual mixing processes and X„ the mass fraction of species 
n. The second term on the right hand side accounts for nuclear 
reactions. The contributions to the diffusion coefficient are con- 
vection, semiconvection, thermohaline mixing, rotationally in- 
duced mixing, and magnetic diffusion. The code includes the 
effect of centrifugal force on the stellar structure, and the trans- 
port of angul ar moment um is also treated as a diffusive process 
(lEndal & Sofi a 1978; P insonneault et al.lll989h . 

The condition for the occurrance of thermohaline mixing is 



5 ^ 



■ Vad < 0, 



(2) 



i.e. the instability operates in regions that are stable against con- 
vection (according to the Ledoux criterion) and where an in- 
version in the mean molecular weight is present. Here (f - 
(dlnp/dlnij)pj, 6 = -(51np/51nr)/>,^, = dlni^/dlnP, 
= (dlnT /din P)ad, and V = dlnT/dlnP. Nume rically 
we treat thermohaline mixin g through a diffusion scheme (iBraunI 
Il997t IWellstein et al.ll200ll) . The corresponding diffusion coef- 
ficient is based on the w ork of [Stern. (.1960) . .Ulrich. (.1972) . and 
iKippenhahn et al.1 (Il980h : it reads 



D„ 



-ath 
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3K 



2 pep (V,d-V) 



(3) 



where p is the density, K - 4acT^ /{3kp) the thermal conduc- 
tivity, and cp = (dq/dT)p the specific heat capacity. The quan- 
tity ath is a efficiency parameter for the thermohaline mixing. 
The value of this parameter depends on the geometry of the 
fingers arising from the instability and is still a matter of de- 
bate jUk ich 1972; Kippenhahn et al. 198C dCharbonnel & ZahnI 
I2007br iDenissenkov & Pinsonneaulti i2008l) . As explained in 
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Fig. 2. Evolution of the region between the H-burning shell 
source and the convective envelope in the RGB phase after the 
onset of thermohaline mixing for a 1 .0 Mq star Green hatched 
regions indicate convection and red cross hatched regions indi- 
cate thermohaline mixing, as displayed in the legend. Blue shad- 
ing shows regions of nuclear energy generation, tracing the H- 
burning shell. 



Sect. l2.2l unless otherwise specified we assume a value ath =2.0 
for the efficiency of thermohaline mixing. This value roughly 
corresponds to the prescription of Kippenhahn et al. (1980), in 
which fluid elements (blobs) travel over length scales compara- 
ble to their diameter. 

For rotational mixing, four different diffusion coefficients are 
calculated for dynamical shear, secular shear, Eddington-Sweet 
circulation and Goldreich-Schubert-Fricke instability. Details on 
the physics of these i nstabilities and the ir implementation in the 
code can be found in iHeger et al.1 ( l2000 h. 

Chemical mixi ng and transp ort of angular mo mentum due 
to magnetic fields (Spruit 2002) is included as in Heg er et al.l 
(2005). The contribution of magnetic fields to the mixing is also 
calculated as a diffusion coefficient (Dmag), which is added to the 
total diffusion coefficient D that enters Eq.[T] 

We compute evolutionary models of 1 .0 Mq, 1 .5 Mq, 2.0 M© 
and 3.0 Mq at solar metallicity (Z=0.02). The initial equatorial 
velociti es of these mo dels were chosen to be 10, 45, 140 and 250 
kms' ( lTassoull2000l) ; we assume the stars to be rigidly rotating 
at the zero-age main sequence. Throughout the evolution of all 
models, the mass-loss rate of .Reimers. (1975.) was used. 



4. Thermohaline mixing on the giant branch 

We compute the stellar models of 1.0, 1.5, 2.0 and 3.0 Mq with 
solar metallicity. The evolutionary calculations presented are the 
same as in Suiis et al. (2008), to which we refer for the details of 
their main sequence evolution. 

The surface composition of low-mass stars is substantially 
changed during the first dredge-up: lithium and carbon abun- 
dances as well as the carbon isotopic ratio decline, ^He and nitro- 
gen abundances increase. After the first dredge-up the H-burning 
shell is advancing while the convective envelope retreats; the 
shell sour ce then enters the chem ically homogeneous part of the 
envelope. lEggleton et al.l (l2006l) and CZ07 have shown that in 
this situation an inversion of the molecular weight is created by 
the reaction ^^He(^He,2p)'*He in the outer wing of the H-burning 
shell in mode ls of 1.0 and 0.9 Mq. This inversion was already 
predicted bv iUlrichI (Il972h . 



4 



M. Cantiello and N. Langer: Thermohaline mixing in evolved low-mass stars 



30 



^ 25 - 



t 20 - 



15 



RGB 



HB 



AGB 



1.0M, 




Luminosity 



1.256x10'° 1.263x10'° 1.270x10'° 1.277x10'° 

time (yrs) 

Fig. 3. Evolution of the surface abundance profiles of the '^C/'^C 
ratio (dotted red line) and ■'He (dashed green line), and of the lu- 
minosity (solid blue line) from the onset of thermohaline mixing 
up to the AGB for a 1 .0 Mq star. 
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Fig. 4. Evolution of the region between the H-burning shell 
source and the convective envelope in the RGB phase after the 
onset of thermohaline mixing for a 2.0 Mq star Green hatched 
regions indicate convection and red cross hatched regions indi- 
cate thermohaline mixing, as displayed in the legend. Blue shad- 
ing shows regions of nuclear energy generation. 



We confirm an inversion in the mean molecular weight in 
the outer wing of the H-burning shell. This inversion occurs af- 
ter the luminosity bump on the red giant branch in the 1.0, 1.5 
and 2.0 Mq models. The size of the ju-inversion depends on the 
local amount of ^He and in the studied mass range decreases 
with increasing initial masfl According to Inequality [2] this in- 
version causes thermohaline mixing in the radiative buffer layer, 
the radiative region between the H-burning shell and the con- 
vective envelope. We emphasize that the extension of the region 
in which the mixing process is active is not chosen arbitrarily, 
but is calculated self-consistently by the code. This is done at 
each time step of the evolutionary calculation by checking which 
grid points fulfill condition|2] This is a major difference between 
models including thermohaline mixing and mo dels where the ex- 
tra mixing is p rovided by magnetic buoyancy ( Busso et al.l2007t 
iNordhaus et a l. 2008; Denissenkov et al. 2009). Indeed for the 
latter a self-consistent implementation is still not available in ID 
stellar evolution codes, and the extension of the extra mixing has 
to be set arbitrarily. 

In our 1 M0 model thermohaline mixing develops at the lu- 
minosity bump and transports chemical species between the H- 
burning shell and the convective envelope (see Fig.|2]i. This re- 
sults in a change of the stellar surface abundances. Fig.[3]shows 
the evolution of the ^He surface-abundance and of the ratio J213, 
qualitatively confirming the result of EDL06 and CZ07, namely 
that thermohaline mixing is depleting ^^He and lowering the ra- 
tio J213 on the giant branch. As already observed by CZ07, the 
surface abundance of '^O is not affected because thermohaline 
mixing does not transport chemical species deep enough the H- 
burning shell. 

Unlike the 1.0 Mq and the 1.5 Mq model, in the 2.0 Mq 
model thermohaline mixing starts but never connects the H- 
burning shell to the convective envelope. This is a direct con- 
sequence of the lower ^He abundance, which results in smaller 
/i-inversion and therefore in a slower thermohaline mixing, ac- 



' During the main sequence of these stars, the pp chain operates par- 
tially burning hydrogen to 'He, but not beyond, into a wide zone outside 
the main energy-producing region. At the end of the core H-buming the 
first dredge-up mixes this 'He with the stellar envelope. Because the 
main sequence lifetime is longer for lower mass stars, these are able to 
produce bigger amounts of 'He. 



cording to Eq. [3] It is surprising that thermohaline mixing, once 
started in the outer wing of the H-burning shell, does not spread 
through the whole radiative buffer layer In fact the H-shell burns 
in a chemically homogeneous region, meaning that no composi- 
tional barrier is expected to stop the instability. The reason is that 
the region unstable to thermohaline mixing moves too slowly in 
the mass coordinates and never catches-up with the quicker re- 
ceding envelope. This situation is shown in Fig. |4l As a result 
no change in the stellar surface composition due to thermohaline 
mixing is observed during the RGB phase of the 2.0 M© model. 

In our 3.0 M© model the H-burning shell never penetrates the 
homogeneous region left by the IDUR Accordingly thermoha- 
line mixing does not occur during this phase. 

In conclusion our models predict that before the He-core 
burning thermohaline mixing is able to change surface abun- 
dances only in stars with M < 1.5 Mq. 

During the RGB evolution the choice of a,h - 2.0, roug hly 
corresponding to the prescription of Kippe nhahn et aTl digS O*) for 
the thermohaline mixing, allows our stellar models to reach he- 
lium ignition without having depleted too much ^He in the enve- 
lope. The presence of leftover ^He allows thermohaline mixing 
to play a role also during a later evolutionary phase, as we show 
below. 



5. Thermohaline mixing on the horizontal branch 
and during the AGB stage 

While CZ07 and EDL07 investigate thermohaline mixing only 
during the RGB, we followed the evolution of our models un- 
til the thermally-pulsing AGB stage (TP-AGB). Indeed a ^- 
inversion is always created if a H-burning shell is active in a 
chemically homogeneous layer, the size of the inversion depend- 
ing on the local abundance of ^^He. This happens not only dur- 
ing the RGB, but also during the horizontal branch (HB) and 
the AGB phase. As a result thermoha line mixing can also oper- 
ate during these evolutionary phases dCantiello & Lange3l2008l : 
IStancliffdl20Tol) . 

Depending on the efficiency of thermohaline mixing during 
the RGB, the ^He can be exhausted at the end of this phase (e.g 
in the models of CZ07). However, stars that avoid extra mixing 
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Fig. 5. Evolution of the internal structure of a 1.0 M© star from 
the onset of thermohaline mixing to the asymptotic giant branch. 
Green hatched regions indicate convection, yellow filled regions 
represent semiconvection and red cross hatched regions indicate 
thermohaline mixing, as displayed in the legend. Blue shading 
shows regions of nuclear energy generation. 



'© 0.6 




E 0.4 - 



1.10x10' 1.13x10' 



1.17x10' 
time (yrs) 



1.25x1 



Fig. 7. Evolution of the internal structure of a 2.0 Mq star from 
the onset of thermohaline mixing to the AGB phase. Green 
hatched regions indicate convection, yellow regions represent 
semiconvection and regions of thermohaline mixing are red 
cross hatched, as is displayed in the legend. Blue shading shows 
regions of nuclear energy generation. 
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Fig. 6. Evolution of the internal structure of a 1.5 M© star from 
the onset of thermohaline mixing to the AGB phase. Green 
hatched regions indicate convection, yellow regions represent 
semiconvection and regions of thermohaline mixing are red 
cross hatched, as is displayed in the legend. Blue shading shows 
regions of nuclear energy generation. 



during the RGB are observed jCharbonnel & Do N ascimentd 
[l998). For these stars the ^He reservoir is intact at He ignition, 
and thermohaline mixing has the potential to play an important 
role during the HB and AGB phases. This is confirmed by the 
evolutionary calculations presented in Sect. l5.1l and l5.2l 

5. 1 . Horizontal branch 

After the core He-flash, helium is burned in the core, while a H- 
burning shell is still active below the convective envelope. In our 
1 Mq model we found that during this phase thermohaline mix- 
ing is present and can spread through the whole radiative buff'er 
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Fig. 8. Profiles of the reciprocal mean molecular weight (l/i-i) 
in the region above the H-burning shell. The plot shows three 
successive times in a 2 M© model during the horizontal branch. 
The black, continuous line represents the model att = 1.13x10^; 
the green, dashed line shows the same model at t = 1.16 x 10^, 
while the blue, dotted line is the profile at t = 1.21 x 10^. 



layer. This is clear in Fig. |5] where thermohaline mixing (red, 
cross hatched region) extends from the H-shell to the convective 
envelope also after ignition of the core He-burning (HB label 
in the plot). Accordingly surface abundances change during this 
phase, as shown in Fig. [3] Here a change of surface abundances 
is also visible after the luminosity peak corresponding to the core 
He-flash. 

Contrary to the 1 Mq model, in our 1.5 and 2.0 M© models 
thermohaline mixing does not change the surface abundances 
during the HB phase. In the 1 .5 Mq model the instability suc- 
ceeds in connecting the H-shell and the convective envelope only 
at the end of the core He-burning (Fig.|6]l, while in the 2.0 Mq 
model this is never achieved (Fig. |7]i. In the latter case ther- 
mohaline diffusion is confined to a tiny layer on top of the H- 
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burning shell, never spreading through the radiative layer (the 
red cross-hatched region in Fig. This is due to a /i-barrier, 
which stops the development of the instability. In Fig. [8] we show 
the profile of 1 /// for the 2.0 model at three successive times 
during core He-burning: the initial peak created by the reaction 
^He(^He,2p)"^He gets smaller, while a dip begins to be visible at 
slightly higher mass coordinate, i.e. at a lower temperature. This 
//-barrier is responsible for stopping the instability; this process 
is discussed in greater detail in Appendix A. 

In the 3.0 Mq model the H-burning shell enters for the first 
time the chemically homogeneous region after igniting He in the 
core. However, in this case also thermohaline mixing does not 
change the surface abundances because is not able to connect 
the H-burning shell with the convective envelope. 
We conclude that in our models thermohaline mixing during 
the HB changes the surface abundances only in stars with M < 
1.5 Mq. 

5.2. Asymptotic giant brancii 

The subsequent evolutionary phase is characterized by two burn- 
ing shells and a degenerate core. The star burns H in a shell and 
the ashes of this process feed an underlying He-burning shell. 
This is referred to as the asymptotic giant branch (AGB) phase. 

During the low-luminosity part of the AGB thermohaline 
mixing works under the same conditions present in the last part 
of the HB phase (see Fig. |5] label AGB). In 1.0 Mq and 1.5 Mq 
models, thermohaline mixing connects the shell source to the en- 
velope. As a consequence surface abundances change, as shown 
for our 1.0 Mq model in Fig. |3] (label AGB). Similarly to the 
RGB and HB phases, no thermohaline mixing is present in mod- 
els with an initial mass higher than 1 .5 Mq. 

During the most luminous part of the AGB the He shell peri- 
odically experiences thermal pulses (TPs); in stars more massive 
than ~ 2 these thermal pulses are associated with a deep pen- 
etration of the convective envelope, the so-called third dredge-up 
(3DUP). In our 1 Mq model we find thermohaline mixing to be 
present also in the TP- AGB. The instability propagates through 
the thin radiative buffer region ("thin" in mass coordinates), and 
reaches the convective envelope. This situation is illustrated in 
Fig. [TT] But there thermohaline mixing only leads to negligible 
changes in the surface abundances. This because of to the very 
short timescale of this evolutionary stage and because most of 
the ^He has already been burned in previous evolutionary phases. 
Overall in our models we found no impact of thermohaline mix- 
ing on the surface abundances of ^He and on the J213 ratio dur- 
ing the TP- AGB phase. On the other hand thermohaline mixing 
can affect the surface abundance of lithium, as we discuss in 

Sect. It] 

We want to stress here that the presence and impact on sur- 
face abundances of thermohaline mixing during the TP-AGB, 
critically depends on the local ^He abundance and on the value 
of the efficiency factor o-th- This is because the local ^He abun- 
dance is related to the previous history of mixing, which in turn 
also depends on the efficiency a^h of the diffusion process. 
We do not know the correct value of o-th in stellar interiors. 
Indeed ffth could also depend on stellar parameters such as ro- 
tation, metallicity or magnetic fields (see Sect. |6]l, and it could 
well be that it changes in the same star through different evolu- 
tionary phases. Therefore our predictions for the changes of sur- 
face abundances due to thermohaline mixing, especially during 
the TP-AGB phase, are strongly affected by these uncertainties. 
Further study is needed to clarify the picture. 



6. Other mixing processes 

6.1. Other mixing processes in our models 

In our 1 .0 Mq and 1 .5 models we found that in the relevant 
layers thermohaline mixing has generally higher diffusion coeffi- 
cients than rotational instabilities and magnetic diffusion. Figure 
|9] clearly shows that rotational and magnetically induced chem- 
ical diffusion is negligible compared to the thermohaline mix- 
ing in our 1.0 M© model. The only rotational instability acting 
on a shorter timescale is the dynamical shear instability, visi- 
ble in Fig. |9] as a spike at the lower boundary of the convec- 
tive envelope. This instability works on the dynamical timescale 
in regions of a star where a high degree of differential rotation 
is present; it sets in if the energy that can be gained from the 
shear flow becomes comparable to the work which has to be 
done against the potential for an adiabatic turn-over of a mass 
element ("eddy") (Heger 1998). However, if present, this insta- 
bility acts only in a very small region (in mass coordinates) at 
the bottom of the convective envelope. As a result thermohaline 
mixing is still setting the timescale for the diffusion of chemical 
species from the convective envelope to the H-burning shell. 

In models of 2.0 and 3.0 Mq thermohaline mixing is less 
efficient due to the lower abundance of ^He. At the same time 
rotational instabilities and magnetic diffusion have bigger diffu- 
sion coefficients, mainly because these models have initial equa- 
torial velocities of 140 and 250 kms"' respectively. Figure [TOl 
shows how during core He-burning rotational mixing and mag- 
netic diffusion become more important than thermohaline mix- 
ing in the 2 model. The radiative buffer layer is dominated by 
the Eddington-Sweet circulation, dynamical shear, and magnetic 
diffusion. Yet the rotational mixing diffusion coefficient is still 
too small to allow the surface abundances to ch ange apprecia- 
bly in this phase, in agreement with results from iPalacios et aP 
(2QQd). The same conclusion is valid for the magnetic diffusion, 
which has the same order of magnitude as the rotational diffu- 
sion in the radiative buffer layer Our models are calculated with 
the Kippenhahn et al., (1980.) prescription for thermohaline mix- 
ing, which implies a smaller diff usion coefficient with respect 
to that proposed bv lUMchl (Il972i) . As a consequence the result 
that thermohaline mixing has in general a higher impact than ro- 
tational mixing and magnetic diffusion in the relevant layers is 
valid regardless of which of the two prescriptions was chosen. 

6.2. Critical model ingredients 

The results described above are obtained with a particular model 
for rotational mixing and angular momentum transport, for 
which several assumptions need to be made. Here we discuss 
the two most important assumptions in the present context. The 
first assumption is that angular momentum transport in convec- 
tion zones can be described with a diffusion approximation and 
a diffusion coefficient derived from the mixing length theory 
(Sect. [3J. The result is near rigid rotation in convection zones. 
Recent 3D hydrodynamic models of rotating red giant con vec- 
tive zones (ISteffen & FrevtagI 120071: iBrun & Pa lacios"2009^ in- 
dicate that the picture may be more complex. [Brun & Palacioi 
(i2009h find radial profiles intermediate between constant angu- 
lar momentum and constant angular velocity, with a large de- 
pendence on the rotation rate of the star. While no general con- 
clusion can easily be drawn from these studies, we may won- 
der how a reduced angular momentum transport efficiency in 
the convective envelope might affect our results. While detailed 
models would be required to exhaustively answer this question. 
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we can expect that a more rapidly rotating base of the convective 
envelope would lead to less shear, and would thus render shear 
mixing in the layers below the envelope less relevant. 

A second crucial assumption is the ado ption of ma: 
netic angular momentum transport according to Spruit ("2002' 
Even though the Spruit-Taylor dynamo has been criticized 
jDenissenkov & Pinsonneaultl 120071; IZahn et aUbOOTh . an effe- 
cient angular momentum transport mechanism Uke the Spruit- 
Taylor dynamo is clearly neede d to understand the observed 
slow rotation of stellar remnants dHeger et al.ll2005l : ISuiis et aTl 
[2008) . While angular momentum transport through grav- 
ity waves has been advocated as an interesting alternative 
(iTalon & Charbonnell2003Ll2008l) . it remains to be demonstrated 
that this mechanism can work to break the rotation of the he- 
lium core in the post-main sequence stages of stellar evolution. 
Therefore, while the reader should be aware of the related un- 
certainties, using the Spruit-Taylor dynamo at this time appears 
reasonable. 



6.3. Interaction of instabilities 

The discussion of the interactions of thermohaline motions with 
the rot ational instabili ties and magnetic fields is complex. In this 
respect ICanutol (1 1999b argues that shear due to differential rota- 
tion decreases the efficiency of ther mohaline mixing. Not only 
iDenissenkov & PinsonneauM (l2008l) claim that rotation-induced 
horizontal turbulent diffusion may suppress thermohaline mix- 
ing. This is because horizontal diffusion (molecular plus turbu- 
lent) may change the mean molecular weight of the fluid ele- 
ment during its motion. They argue that this horizontal diffusion 
is able to halt thermohaline mixing. We think this argument is 
correct in an ideal situation, in which a single blob of material 
is crossing an infinite, parallel slab. Yet in a star the horizontal 
turbulence is acting on a shell, which can be locally approxi- 
mated to a parallel slab with periodic boundary conditions in the 
horizontal direction. This horizontal layer (shell) is rapidly ho- 
mogenized by the horizontal turbulence. Fingers trying to cross 
this horizontal layer are quickly disrupted and mixed. This re- 
sults in a rapid increase of the mean molecular weight fi in the 
shell, so that the region will become unstable to thermohaline 
mixing. A new generation of fingers is therefore expected. But 
the presence of horizontal turbulence is probably making fin- 
gers an unlikely geometrical configuration: blobs that travel a 
small distance before the turbulence is mixing them on a hori- 
zontal layer are more likely. This way thermohaline mixing is 
not stop ped, but only slowed down. This scenario would favor 
the Kip penhahn et al.l (1980) prescription, which actually pre- 
dicts blobs traveling a distance comparable to their size. 

Another intere sting idea has been proposed by 
ICharbonnel & Zahiil (12007 al) . They claim that internal magnetic 
fields can play a stabilizing role, trying to counteract the 
destabilizing effect of the inverse fi gradient. Their conclusion 
is that thermohaline mixing can be inhibited by a magnetic field 
stronger than lO'* - 10^ Gauss. But they warn that their analysis 
ignores both stellar rotation and the spatial variation of B, which 
results in neglecting any possible instability of the magnetic 
field itself (e.g.. Spruit 1999). 

The instability of magnetic fields below the con vective enve- 
lope o f RGB and AGB stars has been discussed bvl Busso et alJ 
(l2007h and |Nordhaus et al. (2008). They argue that dynamo- 
produced buoyant magnetic fields could provide the source of 
extra mixing in these stars. 
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Fig. 9. Diffusion coefficients in the region between the H- 
burning shell and the convective envelope for the 1 .0 Mq model 
during the HB (t= 1.267 x 10'"). The initial equatorial veloc- 
ity of the model is lOkms"'. The black, continuous line shows 
convective and thermohaline mixing diffusion coefficients, the 
green, dashed line is the sum of the diffusion coefficients due to 
rotational instabilities, while the blue, dot-dashed line shows the 
magnitude of the magnetic diffusion coefficient. 
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Fig. 10. Diffusion coefficients in the region between the H burn- 
ing shell and the convective envelope for the 2.0 M© model 
during core He-burning (t= 1.124 x 10^). The initial equato- 
rial velocity of the model is 140 km s"'. The black, continuous 
line shows convective and thermohaline mixing diffusion coef- 
ficients, the green, dashed line is the sum of the diffusion coef- 
ficients due to rotational instabilities while the blue, dot-dashed 
line shows the magnitude of magnetic diffusion coefficient. 



7. Lithium-rich giants 

Lithium is a fragile element, which is destroyed at tempera- 
tures higher than about 3 x lO^K. For this reason it is ex- 
pected that lithium should decrease from its initial value dur- 
ing the evolution of stars. On the other hand, observations 
have shown that about 2% of giants show strong Li lines (e.g., 
Wallerstein & Sneden 1982; Brown et al. 1989). Some of these 
stars even show surface Li-abundances higher than the interstel- 
lar values. 
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For intermediate mass star s a possible solution was proposed 
by ICameron & Fowled (l9TlS, who showed how a net produc- 
tion of 'Li can be achieved during hot-bottom burning (HBB). 
During HBB the convective envelope penetrates into the H-shell 
burning, where 'Be is produced by the pp-chain. In this situation 
the unstable isotope 'Be can be transported to cooler tempera- 
tures by the convective motions, decaying into 'Li in regions of 
the envelope where the temperature is low enough for lithium to 
survive. This results in Li-enrichment at the surface. 

At solar metallic ity stars below ~ 5 Mp do no t experience 
hot-bottom burning dForestini & Charbonnell[T997h . whereas at 
Z = hot-bottom burning is found down to 3 ( Sie ss et all 
For stars avoiding hot-bottom burning, some other mech- 
anism is needed in order to increase the Li surface-abundance. 
A possibility is that some kind of extra mixing connects the H- 
burning shell and the convective envelope, which in the literature 
is often referred to as the cool b ottom process (CBP). The work 
of ICharbonnel & BalachandranI (|2000 ) supports this hypothesis. 
Indeed they found Li-rich stars to be either red giants at the lu- 
minosity bump or early-AGB stars before the second dredge- 
up, in agreement with the idea that some internal mixing oc- 
curs when the H-burning shell enters a homogeneous region. A 
lithium production during the RGB evolution is also supported 
by the recent observations of iGonzalez et alj (l2009h . who find 
lithium enriched red giants at the luminosity bump or at higher 
luminosities. 

lUttenthal er et alj (l2007l) reported the detection of low-mass, 
Li-rich AGB stars in the galactic bulge. Interestingly two of the 
four stars which show surface-Li enhancement present no evi- 
dence for third dredge-up, and thermohaline mixing is advocated 
as a possible source for the extra mixing. 

In our calculations we found that the Li surface-abundance 
is affected by thermohaline mixing during the evolution of low- 
mass stars. While Li is burned during the RGB and HB, ther- 
mohaline mixing has the potential to enhance the Li surface- 
abundance during the TP-AGB phase. To show this, we com- 
puted stellar evolution calculations of the TP-AGB phase in 1 
and 3 Mq with different values of ath- An example of the evolu- 
tion of the Li surface-abundance in the 3 M© model during one 
thermal pulse is shown in Fig.[l2] Our models qualitatively con- 
firm that this instability can enhance the surface Li abundances 
in low-mass AGB stars, even if we can not qua ntitatively repro- 
duce t he high level of enrichment observed by Uttenthaler et alj 
(l2007h . To reach the values of Uttenthaler et al. (2007) a value 
of ath mu ch higher tha n those proposed by [Kippenhahn et al. 
(11980') and' UlrichI (fT97 2) is needed. As discussed in Sect. El a 
quantitative study requires a better knowledge of the efficiency 
parameter for thermoha line mixing gth- 

The observations of lUttenthaler et aP (l2007h show that only 
4 out of 27 galactic bulge stars are Li-enriched. If thermohaline 
mixing is the physical process providing the high Li-enrichment 
observed, we still have to understand why only 15% of the sam- 
ple show this strong enhancement. One possibility is that these 
stars did not experience thermohaline mixing in previous evolu- 
tionary phases. This would leave the ^He reservoir intact, leading 
to a much more efficient mixing during the TP-AGB phase. 
This scenario requires a way to preven t the extra mixing during 
the RGB and HB phases. Charb onnel & ZahnI (l2007ah have pro- 
posed that strong magnetic fields stop thermohaline mixing in 
those red giants stars that are the descendants of Ap stars. They 
call these stars "thermohaline deviant stars". 
Because the fraction of Ap stars relative to A stars (5-10 %), the 
number of red giants that seem to avoid the extra mixing (~ 4%) 
and the observed fraction of Li-enriched AGB stars (15%) are 
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Fig. 11. Evolution of the region between the H-burning shell 
source and the convective envelope during a thermal pulse in a 
1.0 Mq star. Green hatched regions indicate convection, and re- 
gions of thermohaline mixing are red-cross hatched, as displayed 
in the legend. Blue shading shows regions of nuclear energy gen- 
eration. This model is evolved from the zero-age main sequence 
to the TP-AGB with ath =2. 

similar, it may be possible that we are looking at the same group 
of stars at different evolutionary stages. If this is the case, it re- 
mains to be understood why the process that inhibits the mixing 
during the RGB and HB phases is not at work during the AGB. 

A furth er complication arises from the observations of 
iDrake et"al ] (|2002), showing that the incidence of Li-rich gi- 
ants is much higher among fast-rotating objects. They consider 
single-K giants and find that among rapid rotators (ysin/ > 
8kms"') a very large proportion (~ 50%) is Li-rich, in con- 
trast with a very low proportion (~ 2%) of Li-rich stars among 
the much more common slowly rotating giants. Thermohaline 
mixing is not driven by rotational energy, and if any effect 
would be expected, it would be a lower efficiency of t he mix- 
ing with increasing shear and horizonta l turbulence dCanutol 
119991: iDenissenkov & P insonneaulj l2008h . On the other hand, 
an increase in the mixing efficiency with the rotation rate is 
expected if the physic al mechanism behind the extra mixing 
is magnetic buoyancv (IBusso et al.ll2007t iNordhaus et al.]l2008l : 
IDenissenkov et al.]l2009l) . In this case rotation is necessary to 
amplify the magnetic field below the convective envelope. 
Another possibility is that lithium has an external origin, result- 
ing from accretion and ingestio n of planets or a brow n dwarf 
by an expanding red giant (e.g., ISiess & Liviolll99^allbl) . Mass 
transfer or wind accretion in a binary system is also a possible 
scenario. 

The far-IR excess, which is observed in all fast rotating, 
Li enriched giants , is another interesting p iece of the puzzle 
(iDrake et al.ll2002t iReddv & Lamberllliool . While models in 
which some kind of accretion process occurs could explain the 
IR excess, the internal prod uction of lithium can not reproduce 
these observations (but see iPalacios et al.l 1200 ih . We refer to 
iDrake et al.l (l2002h for an accurate review of the proposed mech- 
anism for the formation of Li-rich giants. 

8. Conclusion 

We qualitatively confirm the results of CL07: thermohaline mix- 
ing in low-mass giants is capable of destroying large quantities 
of ■'He, as well as decreasing the ratio J213. Thermohaline mix- 
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Fig. 12. Evolution of Li surface-abundance during one thermal 
pulse in a 3 model. The black, continuous line shows a model 
evolved with ath = 1000; the blue, dotted line refers to the same 
model evolved with a thermohaline mixing efficiency a± - 200. 
In both cases the model experiences third dredge-up. The evolu- 
tion of the star prior to the TP-AGB has been calculated with o-th 
= 2. 



ing indeed starts when the H-burning shell source moves into the 
chemically homogeneous layers established by the first dredge- 
up. At solar metallicity we find that this process is working only 
in stars with a mass below 1.5 Mq. This result is sensitive to the 
choice of the ffth parameter, which regulates the speed of ther- 
mohaline mixing. 

Our models show further that thermohaline mixing remains 
important during core He-burning and can also operate on the 
AGB — including the termally-pulsing AGB stage. Depending 
on the efficiency of the mixing process, this can result in consid- 
erable lithium enrichment. 

Our calculations show that in the relevant layers thermoha- 
line mixing generally has a higher diffusion coefficient than rota- 
tional instabilities and magnetic diffusion. However, we cannot 
address the interaction of thermohaline motions with differential 
rotation and magnetic fields, for which hydrodynamic calcula- 
tions are required. 

In stellar evolution codes thermohaline mixing is imple- 
mented as a diffusive process. This process acts on a thermal 
timescale, but the exact velocity of the motion depends on a pa- 
rameter a±. This parameter is related to the geometry of the fin- 
gers (or blobs) displacing the stellar material and is still a mat- 
ter of debate. The two widely used prescriptions have a param- 
eter Cth that differs by two orders of magnitude. We used the 
[Kippenhahn et al. ( 1980) prescription, even though we also in- 
vestigated the effect of using different values of ath in a few cal- 
culations. .Charbonnel & Zahn C2007b.) used a much more effi- 
cient thermohaline mixing ( UlrichlligT^) . justifying their choice 
on the basis of laboratory experiments of thermohaline mixing 
performed in water, and on the observations of surface abun- 
dances of red giants. 

But the physical conditions inside a star are very different from 
these laboratory experiments, which clearly cannot be used for a 
quantitative study of this hydrodynamic instability. Moreover it 
is not clear if thermohaline mixing is the only physical process 
responsible for the extra mixing, and therefore it is not possible 
to calibrate its efficiency against the observations. 



We argue that is not possible at this stage to firmly identify 
thermohaline mixing as the cause of the observed surface abun- 
dances in low-mass giants (Gratton et al. 2000). In particular the 
long standing ^He problem cannot be considered as solved. 

In agreement with CZ07 we claim that to clarify the picture 
it would be desirable to have realistic hydrodynamic simulations 
of thermohaUne mixing. 
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Appendix A: how to stop thermohaline mixing 

For models of M > 1.5 Mq we found two situations in which 
thermohaline mixing fails to connect the H-burning shell with 
the convective envelope: 

1 . The envelope is receding in mass coordinates and the ther- 
mohaline mixing is not fast enough to catch-up. This situa- 
tion is shown in Fig.|4] 

2. A chain of reactions rising the mean molecular weight can 
create a barrier that stops the mixing process. 

The second scenario occurs depending on the efficiency of the 
reactions rising the molecular weight in the outer wing of the H- 
burning shell. The reactions responsible for creating this compo- 
sitional barrier are ■'He ('*He,y) ^Be (e",v) ^Li ('H,'*He) "^He. 
In the first one ^He and "^He produce ^Be that rapidly decays 
into ^Li. The lithium easily reacts with a proton producing two 
ff-particles; this way the initial molecular weight of 8/9 rises to 
the value 4/3. This occurs in the outer wing of the H-burning 
shell, at a lower temperature with respect to the region where the 
mean molecular weight inversion discussed by EDL06 and CZ07 
is present. As a consequence thermohaline mixing is halted (see 
Figs.|7]and[8]l. 

Therefore two sets of reactions play a role in the evolution of 
thermohaline mixing in the region between the H-burning shell 
and the convective envelope: 



^He-H^He ^ 2iH-i-4He : 1 -» 6/7 



;u: 8/9 ^4/3. 




(A.l) 



(A.2) 



As a consequence, the efficiency in starting and stopping 
thermohaline mixing is regulated by the local abundance of ^He 
and "^He. The first reaction in the chain IA.2I has a lower rate 
than reaction lA.ll at the temperatures in the region of interest. 
However, ■'He (^He,-)/) ^Be depends linearly on the local abun- 
dance of ^He and ^He, while reaction lA. 1 [ depends quadratically 
on the abundance of ^He. Increasing the initial mass of the 
stellar model, th e ratio ^He/"^He decreases in the radiative 
buffer layer (e.g.. lBoothrovd & Sackmana.1999,) . Consequently 
if one increases the initial mass, the set of reactions IA.2I at 
some point becomes more important than lA.ll meaning that a 
compositional barrier is created. This explains why we observe 
a threshold in mass at about 1 .5 Mq, above which thermohaline 
mixing is not efficient anymore. 

Note that both scenarios 1. and 2. depend on o-th, i e- on the 
speed of the mixing process. Increasing the value of ffth makes 
thermohaline mixing more difficult to be stopped; therefore the 
value of the mass threshold depends on the choice of ath. 
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